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Abstract 

It has been shown recently 10 that Cauchy transforms of orthogonal polynomials appear nat- 
urally in general correlation functions containing ratios of characteristic polynomials of random 
N X N Hermitian matrices. Our main goal is to investigate the issue of universality of large N 
asymptotics for those Cauchy transforms for a wide class of weight functions. Our analysis covers 
three different scaling regimes: the "hard edge", the "bulk" and the "soft edge" of the spectrum, 
thus extending the earlier results known for the bulk. The principal tool is to show that for finite 
matrix size N the auxiliary " wave functions" associated with the Cauchy transforms obey the same 
second order differential equation as those associated with the orthogonal polynomials themselves. 
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1 Introduction 



Characteristic polynomials are important objects in Random Matrix Theory (RMT) and enjoy a wide 
range of applications in various fields of physics and mathematics. To name only a few, we mention 
the applications to Quantum Chromodynamics (QCD) jl], Quantum Chaos |2| as well as their role in 
elucidating properties of the Riemann zeta function along the critical line [Hill]. In fact, products of 
ratios of the characteristic polynomials at different values of the spectral parameter provide the most 
general generating functional for all spectral correlation functions in RMT. Characteristic polynomials 
are also useful in providing a direct link between the RMT description and the effective field theory 
such as Chiral Perturbation Theory (xPT) in the application to QCD [Slinil3- There, products of 
the ratios of characteristic polynomials enjoy a direct interpretation as partitions functions including 
fermionic and bosonic degrees of freedom. However, these objects are very difficult to obtain in full 
generality from the finite volume limit of xPT, where one has to parameterize and calculate supergroup 
integrals jS]. Alternatively, one has to perform replica calculations exactly |7j, a possibility opened 
due to recent progress in understanding the analytical structure of the moments [H]. From this side, 
correlation functions of characteristic polynomials recently calculated within RMT [^1 ^] provide the 
result yet to be obtained [llj directly from the field-theoretic calculations. 

An important issue for all RMT calculations is the question of universality. By this one usually 
understands that the predictions of RMT for a given object are insensitive to the details of the 
particular models. When this property holds, it must be a consequence of the common underlying 
symmetry of that class of the models. This type of argument validates applying RMT results to 
systems of varying nature in all the diverse situations mentioned above. 

The character of universality depends essentially on the scale in the spectrum one is interested in, 
as well as on the nature of the singularity at the support of the spectrum (the so-called "edges" of 
the spectrum). In the following we will be only concerned with the scales comparable with the mean 
level spacing. This issue is widely known as that of "microscopic" universality. Suppose first one is 
interested in such quantities as the correlation function of eigenvalue densities, or the distribution of 
the largest/smallest eigenvalue, or some related objects. Then restricting oneself to the class of random 
matrix measures invariant with respect to the change of basis (here the unitary invariant ensembles) 
the quantity of interest can be expressed as a determinant containing kernels made of polynomials 
orthogonal with respect to that measure In this way studying universality amounts mainly to 
investigating asymptotics of the orthogonal polynomials in various regimes, and is well established by 
now both on physical [HI [HI HSl [El [13 CHI IH as well as mathematical [SOI EH |22l EHl level of rigor. 

On the other hand, more general correlation functions containing products and ratios of the char- 
acteristic polynomials (in particular, those used to generate all the eigenvalue correlations by differen- 
tiation and taking the discontinuity of the resolvents) attracted considerable attention only recently 
[3 ini ^1 ^3 [211 • While for the products only the universality follows quite straightforwardly from 
the mentioned results [^ El 1^ [101 [^ HH HB] (see also |23j ) , the universality of expressions 
containing inverse powers or ratios was properly understood only in the bulk of the spectrum f25'. 

There are several possibilities to calculate correlation functions containing ratios of characteristic 
polynomials. For the case of the Gaussian measure the only tool until quite recently was the exploita- 
tion of anticommuting Grassmann variables (the method of " supersymmetry" due to Efetov or its 
variants .27^ 28.. ^ [HI El ) • An alternative for general unitary- invariant measures was discovered quite 
recently by Fyodorov and Strahov |in[ I25j and resorts again to the orthogonal polynomial method. 
An important new feature, however, is that as soon as inverse powers of characteristic polynomials are 
taken into consideration, one has to deal not only with the orthogonal polynomials but also with their 
so-called Cauchy transforms. More precisely, a very general result - valid for all invariant ensembles of 
finite-size complex Hermitian matrices |10| - reduces arbitrary correlation function containing ratios of 
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characteristic polynomials to a determinant made of orthogonal polynomials and their Cauchy trans- 
forms. The universality of the corresponding expression in the microscopic limit in the "bulk" of the 
spectrum was rigorously demonstrated by Strahov and Fyodorov 25, by using the Riemann-Hilbert 
approach. Those authors however considered neither the "soft" nor the "hard" edge of the spectrum, 
the latter being a characteristic feature of the so-called " chiral" version of the unitary ensemble |291 ^ . 

Our goal in the present paper is to employ an alternative method and to obtain universal asymp- 
totics for all three regimes, including spectral edges of both types. This is particularly important in 
view of the application to QCD and the role of characteristic polynomials played there (see. e.g. in 
[7| lllj and references therein) . The method we are going to exploit to extract universal features of the 
correlation functions is based on the works by Kanzieper and Freilikher [HHl- The clear advantage of 
this technique is its explicit form and conceptual simplicity, although it lacks the rigor and vigor of the 
Riemann-Hilbert technique (and in fact makes use of some facts proved by the latter). First we will 
show the existence of an exact second order differential equation satisfied by the Cauchy transform 
of the orthogonal polynomials for a large class of measures. Then by considering the microscopic 
large- limit of that differential equation we can demonstrate their universality at the origin ("hard 
edge"), in the bulk and at the soft edges of the spectrum. The procedure follows exactly the same 
lines as for the orthogonal polynomials themselves (see j3()|). since in both cases we have to deal with 
the very same equation. The universality of asymptotics for the general expressions ^Uj of arbitrary 
ratios of characteristic polynomials then follow. It is necessary to mention that the method implicitly 
uses an assumption of the existence of a single- valued function arising in a certain limit for recursion 
coefficients relating orthogonal polynomials of different order. This was indeed shown to be true in the 
case of a spectrum supported on a single interval in the framework of the Riemann-Hilbert technique. 

The article is organized as follows. In sect. |2lwe give the main definitions and recall the main result 
of ^Uj. In the following sect. |21we derive a differential equation for the wave functions associated 
with the Cauchy transforms for finite and show that it is identical to that for the wave functions 
of orthogonal polynomials. For convenience we treat the unitary ensemble (UE) first and translate to 
the chiral unitary ensemble (chUE) later on. Sect. 0] deals with the large- limit yielding universal 
results for both unitary and chiral unitary ensemble in all three microscopic regimes. We close with 
the discussion and open questions in sect. [3 

2 Definitions 

In this section we introduce the main objects we wish to study in the present paper. The matrix model 
partition function with respect to which we calculate expectation values is given by the following 
integral over the (unitary invariant) Haar measure for matrices H: 



The matrix H can be either complex Hermitian, representing a general unitary invariant ensemble 
without chiral symmetry, or it can be decomposed into a product of complex rectangular matrices 
W of size N X (N + u), H = W'^W, representing the chUE. In the first case the matrices H can be 
diagonalized by a unitary transformation with real eigenvalues A = diag(Ai, . . . , Aat) are distributed 
according to a weight function w{X) with support on M. In the second case we deal rather with real 
positive singular values of the matrices W with the support being M+, respectively. We have also 
employed the notation A7v(A) = ni>j(-^« ~ -^i) the Vandermonde determinant. While the results 
cited in this section hold for general weight functions with convergent moments on some domain 
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Z? C M, we will concentrate mainly on weights of the form w{\) = exp[— A^T^(A)] in the following 
sections. The potential V{x) may contain both an (even) polynomial part and a logarithmic part. 

We will be interested in extracting the large-A^ universal behavior of the correlation functions 
defined via expectation values of arbitrary ratios of characteristic polynomials: 



}CN{Mf,Mb) 



L 



n det (mjljv + -H' 



M 



(2.2) 



In the context of QCD the matrices Mf = diag 



m 



f 



Mh = diag {m\, m^) are said to 



contain the fermionic (/) and bosonic (6) mass terms. The expectation value of ratios of characteristic 
polynomials is known as the supersymmetric partition function. It may be used to calculate correla- 
tion functions of the eigenvalue densities and serves as a general generating function for multi-point 
resolvent operators. 

In order to give rise to positive definite partition functions in the case of the UE the masses have 
to be chosen purely imaginary and occur in pairs of complex conjugates. In this way only positive 
definite terms appear, det(i?Ti + H) det{—im + H) = Ylki^'^ + ^k)- chUE each determinant 

factor is positive provided we parameterize det(m^ + W^^M^) = Hfe^C"^^ + ^k), > 0. At the same 
time, from the mathematical point of view there is no need to impose any restrictions on masses in 
eq. (|2.2() . apart from requiring infinitesimal imaginary parts for all m^. 

As it has been shown in ^01 for general weight functions and finite size the generating function 
eq. 1)2. 2|) can be written as a determinant of size M + L containing polynomials orthogonal with 
respect to the weight w{X). Let us define the monic orthogonal polynomials 7r„(A) = A" + 0{X^~^) 



ciAu;(A)7r„(A)7rm(A) = 5„m/n 



(2.3) 



D 



with normalization factors /„ uniquely defined through the Gram-Schmidt orthogonalization.^ The 
second building block are the Cauchy transforms of those polynomials which appear only as long as 
M > 0, 

Me) ^ ^ I d\ ^7r„(A) , 9m(6) / . (2.4) 
The result derived in ^0] then reads ^ 



}C{Mf,M,) 



oc 



A{Mh)A{Mf) 



det 



'&N+L~i {-m\) 



"Sn-m {-nAi) ■&N-M+1 {-nAj) . . . '&N+L-1 



ttn-m 



TTN-M 



-m. 



T^N-M+l 



T^N-M+l 



-m. 



T^N+L-l 



-mi 



-m, 



. (2.5) 



^For weight functions on D = R the corresponding norms can be calculated recursively solving nfn-\ = 
— J^^dX w'(A)7r„(A)7r„_i(A) for even potentials plus an additional equation for potentials without parity. 

■^The object eq. 12.51 can be given a meaning for real values of the mass parameters mj" as well, by choosing the 
principle value in eq. I|2.4^ . However, to use it as a generating functional we are precisely interested in its discontinuities 
in the complex plane. 
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3 A differential equation for the orthogonal polynomials and their 
Cauchy transforms 

As is well known I3UI 02] it is convenient to relate to the system of orthogonal polynomials vrn(A) 
a system of associated objects: 



tpnW = exp 



N 

-yK.(A) 



(3.1) 



In 



frequently called the "wave functions". The name has its origin in two properties: (i) the functions 
ipniX) obey a closed second order differential equation for finite which can be interpreted as a kind 
of one-dimensional Schroedinger equation and (ii) the functions are orthonormal. The first property 
can be shown, in particular, by the formalism introduced by Kanzieper and Freilikher [20], to which 
we refer an interested reader for a detailed review. Below we exploit the same formalism to infer 
properties of the Cauchy transforms of the orthonormal polynomials, after suitable redefinition of 
the wavefunctions. This will enable us to treat the Cauchy transforms on equal footing with the 
polynomials themselves when extracting the universal properties in the microscopic large-A^ limit. 
The wavefunctions of the orthogonal polynomials and of their Cauchy transforms will be then two 
independent solutions of the second order differential equation, with the latter being singular at the 
origin. 

Due to an intimate relation between the polynomials of the UE and the chUE for general weights 
w{X) we can restrict ourselves to treat only the UE in greater detail in the following subsection. Once 
our results are established those for the chUE will follow rather straightforwardly, as discussed in in 
the subsection 13.21 

3.1 The Unitary Ensemble 

For the UE the eigenvalues of the complex Hermitian matrix H are distributed along the full real line. 
To be specific let us introduce the following parameterization of the weight function 

w{\) = A2"exp[-W(A)] = exp[-NVa{X)] ■ (3.2) 

Here we have introduced the potential 

2 '^^^ 

V^{X) ^ V{X) - ^ ln(A) , V{X) = ^ ^-^X'' , (3.3) 

k=l 

with an even polynomial part V{X) of arbitrary but fixed even degree d. In order to apply the 
formalism of [30] we switch from monic polynomials 7r„(a;) to orthonormal polynomials. 



6m = / dXw{X)pk{X)pi{X) , pk{X) = T:k{X)Hfk ■ (3.4) 
The potential in eq. (|3.3|) being even, the polynomials possess the following parity property: 

Pn{-X) = i-ir Pn{X) . (3.5) 

As is well known (see e.g. j221 l3Uj and references therein), all orthogonal polynomials on the real line 
or a subset of it satisfy a three-step recursion relation 

Xpn{X) = Cn+lPn+l{X) + (A) , (3.6) 
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where the middle term vanishes because of the parity. The recursion coefficients can be obtained from 
the normahzation factors /„ of the polynomials as: 



(3.7) 



y Jn-l 

The second operation, which together with the multiplication by A in eq. (|3.6|) closes on the ring of 
orthogonal polynomials, is differentiation over A. Following it can be most conveniently written 
as 

K(A) = An{X)pn-l{X) - Bn{\)PnW ■ (3.8) 

Here, the functions An{\) and Bn{\) are defined as 

An{\) = N Cn r dll^K{n)-VLW]Pn{ri? 
J-oo - A 



Ncn dii ^[^'(r?) - V'{\)] Pn{riY , (3-9) 
r] — X 



— oo 



J-oo V-X 

^ Bi:^{X) + (l-(-l)")^. (3.10) 

Because of parity the singular a-dependent term drops out in ^^(A). Inserting the explicit form of 
the potential from eq. 1)3. 3|) it can be easily seen that A„(A) and the regular part Bn^{X) of Bn{X) 
are polynomials of degree d — 2 and d — 3 respectively. Prom their definition one can easily show the 
following identity to be valid jl5j : 

S„(A) + S„_i(A) + iVV;:(A) = -^An-i{X) . (3.11) 

Cn-l 

This identity, when combined with the two properties of the polynomials - the recursion relation eq. 
()3.6() and the differentiation eq. 1)3. 8|) - suffices to show that the following second order differential 
equation is satisfied by the wavefunctions defined in Eo. p.ll) :. 



where 



<(A)-F„(A)V';(A) + G„(A)^4A) = 0, (3.12) 



F„(A) . (3.13) 



G„(A) 



^^„(A)^„_i(A) - fi?„(A) + ^Ki(A)) 

+ {b^W + !^KiX))' - ^ (b..(a) + !iv:.(X)) . (3.14) 



Here ' denotes the derivative with respect to A everywhere. The differential equation Ea. H3.12() holds 
for finite-A^ and arbitrary weight w{X) corresponding to a given polynomial ^(A) of degree d and 
parameter a > —1/2. The fact that a can be real will enable us to translate most easily the results 
for the UE obtained below to the chUE in the next subsection. 
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Our goal is to derive a similar differential equation for the Cauchy transform at finite- A'^. In terms 
of the orthonormal polynomials they are given by 



+~ w(X) 

dX Y^Pn(A) , 9m(e)/0 . 
00 A — e 



(3.15) 



This defines in fact two functions, one in the upper and one in the lower half plane for positive or 
negative imaginary part, respectively. Each of it can be extended to the other half plane using the 
relation 

hn{-e) = (-1)"+! K{e) . (3.16) 

which follows easily from the parity of the polynomials, eq. ()3.5() . We note that the Cauchy transform 
has the opposite parity than the corresponding polynomial. Let us stress that the transformation 
()3.16l) does not relate the two different functions defined in eq. ()3.15|1 . 
In the following we show that the object 

K{e) = exp[+NVai€)] K{e) (3.17) 

satisfies the same recursion relation eq. (|3.()|) and differentiation relation eq. (|3.8|) as the corresponding 
orthogonal polynomial Pn(A), provided that n > d—1. As an immediate consequence the wavefunction 
for the Cauchy transform defined as 



Xn{e) 



exp 



N 



exp 



N 



(3.18) 



satisfies the second order differential equation identical to Eq. 1)3.12(1 at finite-n, provided the condition 
n > d — 2. The corresponding solution of the differential equation at finite and large-n will then have 
to be chosen according to the appropriate boundary condition for ipn and Xn, respectively. 
Let us first consider multiplication on /i„(e) for n > 0: 

1 



e/i„(e) 



2TTi 

-1 
2tH 



dX- epn{X) 



00 
+00 



A-e 
dX w{X)pn{X) + 



1 
27d 



+°° w(X) 

dX- [Cn+lPn+l{X) + (A)] 

■00 A — " 



c„+i/i„+i(e) + c„/i„_i(e) , n > . 



(3.19) 



At the first step we have added and subtracted A to cancel the denominator and used the recursion 
relation eq. (|3.6|) . At the second step we used that Pn(A) is orthogonal to po{X) = const for n > 0. The 
fact that the recursion Eq. (|3.19|) also holds for hn{e) as well as for Xn(e) is an immediate consequence. 
The next step is to consider differentiation acting on /in(e) for n > d — 2: 



deKie) 



1 

1 

27ri 



+00 



dX 



,N{V^ie)-V^iX)) 



—00 
+00 



A 



[N{Vi{e)-ViiX))pn{X) + K(A)] 



dX QN{Vo,{e)-Vo,{X)) 



N 



A-e 



Pn(A) - ^Pn(A) 



eA^ 



^^n(A) i3-g(A) + (l-(-l)")«A-i ■ 

+ ^ Pn-l{X) r Pn{X) 



X 



+00 



— / dX e^(^-W-^-W) 

2-711 



Blt^ie) + (1 



-l)")ae-i 



A 



A 



2a Anie) 
rPn(A) + -^p„_i(A) 

£A A — £ 

a(i-(-i)") 

Pn(A) + ^ ^^Pn(A) 



eA 



Anie)K-iie)- i3-g(e) + (l-(-l)")- /i„(e) , n>d-2 



(3.20) 
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Here we first have used the relation dejr\ = —dx-^:rx performed the integration by parts. The 
second step spells out the definition of Va{X) and exploits Eq. H3.8|) for p^(A). We can now use that 
{V'{e) — V'{X))/{X — e) is a polynomial of order d — 2 and thus orthogonal to Pn(A) for n > d — 2. 
The same argument helps to rewrite ^„(A)/(A — e) = (^„(A) — A„(e))/(A — e) + A„(e)/(A — e) and to 
drop the first term, which is again a polynomial (this times of the order d — 3) and thus orthogonal to 
Pn{X). We can repeat the same steps for Bl^^(\) and the poles 1/A. The remaining terms ~ l/(Ae) 
require a little bit more careful treatment. For n even pn{X)/X is odd and thus the corresponding 
integral vanishes identically. For n odd the last term (1 — (—1)") gives a contribution that precisely 
cancels the first one. We thus finally arrive at the same relation for differentiation of hn{e) as that 
given by eq. p?8|) . Consequently, Eqs. (|3T9|) . (|3?20|) and p?T8|) imply 

X':,{X) - Fn{X)x'n{X) + Gn{X)Xn{X) = . (3.21) 

Thus, the wave functions associated with the orthogonal polynomials and their Cauchy transforms 
satisfy the same differential equation, provided that condition n > d — 2 holds. This constraint is not 
an obstacle for large but finite-n, for any given potential V{X) with a fixed finite degree d. 

The virtue of our analysis at finite n is that it obviously translates all the universality proofs 
carried out at the origin, in the bulk and at the edge of the spectrum for the orthogonal polynomials 
to their Cauchy transforms. We will concentrate on the the large-n limit in the next section 01 but 
first we repeat the above analysis for the chiral case of chUE. 

3.2 The chiral Unitary Ensemble 

The eigenvalues for chUE are distributed on the positive real line and governed by the following 
chiral weight function: 

w^''{X) = A"exp[-iVT/'='^(A)] = exp[-W„"'^(A)]. (3.22) 
Here the potential V^^{X) is of degree d/2, with d even, and is defined as 

d/2 

V^\X) ^ ^MA'=-|ln(A). (3.23) 

k=l 

The orthogonality relation for the corresponding orthogonal polynomials denoted as P^'^hW reads 

/•oo 

6ki = / dX w^\X)pi%iX)pl^^,{X) . (3.24) 
Jo 

Here we explicitly display the dependence on the parameter a. Substituting X^ = y one can easily 
see that the orthogonal polynomials of the chUE can be obtained from those of the UE by simple 
replacement of the argument [18j : 

Pl%i^') = piT'^^'W ■ (3.25) 
Here we have to appropriately identify the corresponding potential as 

V^Hy') = 2 V{y) , (3.26) 

which effectively means that the non-chiral polynomials are taken at values 2N. In fact, one can look 
at the relation Eq. (|3.25|) as a generalization of the well-known relation between generalized Laguerre 
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and Hermite polynomials. Consequently, there is no need to repeat the whole analysis of the previous 
subsection but one can simple read off the chiral orthogonal polynomials and their asymptotics from 



the same differential equation (|3.12)) for even at shifted value of a. 

Let us now show that the Cauchy transforms of the chiral orthogonal polynomials relate to those 
of the non-chiral ones as well. We first define the Cauchy transform for the chiral case as 

Substituting \^ = y and inserting the relation ()3.25|) we obtain 

= ^/_r^^^«"'^^"'^^^^^i':^c.(^^) (3.28) 

where we have used the parity of the non-chiral polynomial P2^^^'^\y)- On the other hand we can 
rewrite the Cauchy transform of the non-chiral ensemble eq. 1)2.4(1 with an even index 2n as 

4? (^) = ^ dy-^^ e-^^y^^y) pt\y) . (3.30) 
Comparing eqs. ()3.29|) and ()3.30j) together with eq. ()3.26|) we thus arrive at 

hi^U^') = 7 ht^'^'\e) , (3.31) 

where the right hand side has to be taken at 2N as for the polynomials themselves. Consequently the 
same analysis for the "chiral" Cauchy transforms follows from the same differential equation ([3.21(1 . 

We note that while the Cauchy transform eq. (|3.27|) defines a single function with a cut along the 
positive real line M+, the second transform in eq. ()3.28|) defines again two functions in the upper and 
lower half plane as in eq. ()2.4|) . It is in this latter picture on R that the results of the UE can be 
translated to the chUE. This is why in the following section we will always give two different solutions 
for the chiral Cauchy transforms as well. 



4 Universality 

In this section we demonstrate the universality of the asymptotic polynomials and their Cauchy 
transform in the microscopic limit, which then translates to all correlation functions of characteristic 
polynomials through eq. ()2.5|) . The limit is called microscopic since it includes a certain rescaling 
of the eigenvalues with an appropriate power of the parameter N, according to the spectral region 
we investigate. We will first study the microscopic limit at the origin (frequently called the "hard 
edge"). This case is singled out due to the presence of logarithms in the potential. The correlations 
in the bulk of the spectrum then follow easily. Furthermore we also study the edge of the spectrum 
defined by the endpoints of the support of the macroscopic or smoothed spectral density at large N 
("soft edge"). Here we will restrict ourselves to the cases where the spectrum support is a simply 
connected interval, the so-called one-cut phase. We will also stay away from multicritical points, 
where the macroscopic density develops extra zeros at the boundary of spectra or inside the support. 
This amounts to requiring the limiting functions for A]\[{X) to be nonvanishing on the support as well 
(see eq. ()4.3j) below). While such multicritical cases can be, in principle, analyzed in the framework 
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of the same approach ^1 El 1^ the issue becomes much more subtle due to subleading terms in the 
large- A'^ expansion. 

Our first step is to derive the large- limit of the differential equations (|8.12j) and (|3.2H) . Our 
assumption of the one-cut support is equivalent to the fact that in the large- limit the recursion 
coefficients c„ in eq. p.Tf) approach a single valued function 

lim Cn = c{t = n/N) . (4.1) 

While in the earlier physical matrix model literature this assumption was introduced as an ansatz 
supported by numerical evidence, it has been proved rigorously in the recent mathematical literature 
PT] . Eq. (|4.1() implies the following asymptotic form for the functions ^Ar(A) and -Bat (A) in the 
differential equation pUj : 

= 7rp(A)^=L= , \€ {-a,a) . (4.3) 

Here a = limjv^oo 2cAr is the endpoint of the spectral support. For convenience we have also presented 
the relation to the macroscopic density of eigenvalues denoted by p{X). Using the recursion relation 
eq. H3.6|l it can be easily seen that ^(A) remains a polynomial of degree d — 2, but the coefficients get 
smoothened by taking the large- limit. This corresponds to the one cut phase through eq. (|4.3() . 
with the nonpolynomial part of p(A) ~ — A^ being a simple square root. The regular part Bn^{X) 
of the function Bn{X) also remains in the same limit a polynomial of order d — 3, with the limiting 
form jl5j : 



(4.4) 



Using these expressions we obtain the following result for the combination entering the differential 
equation eq. (|3.11|1 

l\m^ (^B^{X) + y^^(A)) = AT ^A{X) - (-1)^ ^ . (4.5) 

We have kept the last term as it will give rise to two different equations for even and odd A^. Con- 
sequently, Gat (A) can be expressed asymptotically only in terms of A{X) and therefore via p(A) by 
virtue of Eq. ()4.3() . We have: 



Although it seems that the first term in Gat (A) is always dominating the second term may contribute 
in some regimes as well, depending on the chosen rescaling of the eigenvalues which we are about to 
investigate in the next subsection. 
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4.1 The hard edge 

We start by defining the microscopic rescaling in the large- Hmit close to the hard-edge origin of the 
spectrum to be 

NX = ^ . (4.8) 

This amounts to magnifying the correlations at the point A « as ^ is kept fixed while — > oo. 
Inserting this rescaling into the differential equations H3.12() and H3.21() and using the asymptotic 
expressions for the coefficients, eqs. (lOll and (jOI) it is easy to see that the two terms with the 
logarithmic derivative of A(X) are both of subleading order 0{1/N). We thus obtain 

V^NiO + (vrV(O)' + (zil^IjH^^ ^^(^) = , v'^v G {i^N, xn} • (4.9) 

This expression is clearly universal as the dependence on all the specific parameters {g2k} in the 
polynomial potential in eq. (||-{.2j) has been translated into a single parameter - the macroscopic spectral 
density at the origin p{0). Due to the oscillating sign we have to distinguish between even and odd 
values of N, corresponding to even and odd polynomials and their Cauchy transforms. The second 
order differential equation has two independent solutions and it is the requirement of regularity of the 
solution at the origin that dictates which one to be chosen. The wave functions associated with the 
polynomials iPn{0 regular at the origin, fixing the corresponding solution to be Bessel functions of 
the first kind Ju{£,)- the other hand the wave functions associated with the Cauchy transforms are 
singular at the origin. As we will see below this makes us to choose solutions to be Bessel functions of 
the third kind also known as the Hankel functions, h!)'\Q or H^\q, depending on the sign of the 
imaginary part of the complex argument A^e = We therefore obtain 

V'2n(0 ~ (7rp(0)O^J,„i(vrp(0)e) , 
^2n+i(0 ~ (vrp(0)e)^J,+ i(vrp(0)e) , (4.10) 

(7rp(0)C)Wi'\(vrp(0)C) , 9m(C)>0 



""^^^ ^ -(..(0)C)^i/f :(vrp(0)C) 



3=m(C) < 



(7rp(0)C)^i/^';i(vrp(0)C) , 9m(C)>0 
"^"-^^^^^ ^ ^ -(MO)C)WS:(vrp(0)C) , 3m(C)<0 " ^''''^ 

In order to extract from these wave-function asymptotics the expressions for the bare polynomials 
and their Cauchy transforms we still have to remove the corresponding weight factor from eqs. (|3.1j) 
and H3.18() accordingly. The microscopic limit of the weight function eq. 1)3. 2() reads: 

lim N'^'^w{X) = . (4.12) 

Multiplying or dividing the respective wavefunction by the square root of the weight function, we 
arrive at the final universal expressions to be inserted into the determinant of eq. 1)2. 5(1 : 



P2n(6 ~ r"+^^,_i(vr/5(0)0 

2 



P2n+m ~ r"+5j ,(vr/,(0)O , (4.13) 
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/i2n(C) 



2n+l 



(C) 



r+5i/«,(vrp(0)C) 

" 2 

r+ii7«,(^p(o)c) 



3=m(C) > 
9m(C) < 

9m(C) > 
9m(C) < 



(4.14) 



Let us now justify our choice of the type of Bessel functions featuring in eq. (|4.11j) . Due to the 
singularity at the origin we have to choose between Bessel functions of second kind (also known as 
Weber functions), i^(C)) and the Hankel functions. We give two independent arguments in favor of 
the latter choice. The first reason for selecting the Hankel functions is that it allows one to match 
the known "bulk" asymptotic result |25j after taking the asymptotic limit ^ ^ oo. While the Weber 
functions would asymptotically produce sine and cosine terms (as the Bessel functions of first kind 
do), the Cauchy transforms in the bulk were shown to approach exp[ibiC] jj^]. This asymptotic is 
correctly reproduced by the Hankel functions. The second argument is related with the necessity to 



distinguish between the two types of the Hankel functions h!)^^ {(^) and Hl"^' according to the sign 
of 9m(C). 

In fact, for a = and Gaussian potential w{X) = exp[— A^A^] one can explicitly write down 
the integral representations for the polynomials and Cauchy transforms at finite-A^, as discussed 
in Switching back to the monic normalization the polynomials are the Hermite polynomials 

7r„(A) = iJ„(A\/iV)/(2^/iV)". It follows that 



r(2)< 



7r2m(A) 
Vr2m+l(A) 



POO 

-ir2NVX ds e-^''s'^'^+2j_,(2NsX) 
Jo 2 

POO 

-ir2NVX dse-^'%'^"'+^J^,(2Ns\), (4.15) 
Jo 2 

where we have used the standard integral representation of the Hermite polynomials as well as re- 
expressed cosine and sine by the half-integer Bessel functions J_i_i. As a consistency check one can 
verify that eq. 1)4. 13() at a = follows from a saddle point approximation in the microscopic large- A^ 
limit, remembering that vr/9(0) = \/2 for the Gaussian potential. Using eq. (|4.15|) we obtain a finite- A^ 
integral representation of the Cauchy transforms as follows. We can rewrite 

1 f°° 

= sgn[9m(e)] i dt ex.p[—sgn[Qm{e)]it{x — e)] , 3'm(e) / . (4.16) 

x-e Jo 



Inserting eqs. (|4.15() and 1)4. 16() into the definition ()2.4|) we can first perform the integral over A and 
then the one over t, arriving at 



1? 



2m+l 



H^y,{2Nse) 
-H^^l{2Nse) 

2 

H^^l{2Nse) 
-H^\{2Ns€) 



9m(e) > 
9m(e) < 

9m(e) > 
9m(e) < 



(4.17) 



In taking the microscopic large-A^ limit we again recover eq. ()4.14|) at a = from a saddle point 
approximation. This fact confirms the correct choice of the solution for the Cauchy transform and 
thus for the wavefunctions eq. 1)4.11(1 for a Gaussian potential with a = 0. Since the parameters a as 
well as the higher order coupling constants g2k in the potential eq. ()3.3|) are real and can be switched 
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on smoothly the choice of the type of Hankel function remains fixed for a general potential for general 
a as well. 



Now we can use the relations between the polynomials and Cauchy transforms of the chiral and non- 
chiral ensemble, eqs. (|3.25() and (|3.31|) respectively, to read off immediately the universal asymptotic 
formulas for the chiral Unitary Ensemble from H4.13j) and H4.14p : 

v^ScH^i^) - r°^(^p(o)e), (4.18) 

..2^ / C"^i'^(vrp(0)C) , am(C)>0 



We see that for integer values of a the indices of the Bessel functions are also integer, compared to 
half integer values for the non-chiral case above. To fix independently the correct Bessel solution 
corresponding to the sign of 9m(^) we can repeat the argument as in eqs. (|4.15() and ()4.17|1 . but 
this time for general a 7^ 0. The corresponding finite- expression have already been given in |llj . 
following the integral representation of generalized Laguerre polynomials, and we only repeat them 
here: 

^...(e^) = t^Ne'^ r dt.--H---4 ^Sr^^l ' ^"^^'^^^ (4-20) 
' 2 Jo \ -H^^^ {2NeVt) , 9m(e) < ^ ' 

The saddle point approximation reconfirms eq. (|4.19jl for the Gaussian potential V^^{X) = A — 
and, by invoking the continuity argument for a general potential as well. 



4.2 The bulk of the spectrum 

Li this section we briefly describe the microscopic limit in the bulk of the spectrum, which amounts 
to setting a = in the differential equation ()4.9() . For the non-chiral case the corresponding formulae 
match with those obtained in |25j . 

It is evident that when logarithmic terms are absent the origin of the spectrum is no longer singled 
out. Under these conditions correlation functions when properly scaled should be the same at every 
point inside the support of the spectrum, due to translational invariance. More precisely, let us pick 
a fixed value A^, € {—a, a), that means neither at the origin nor at the edge of the support. The 
microscopic rescaling around the chosen point is then defined as 

iV(A-A,) = ^ . (4.21) 

Inserting this expression into that for the asymptotic coefficients, eqs. ()4.6p and 1)4. 7() . we can easily 
see that only the first term in Gat (A) contributes yielding 

</''^(e) + vrV(A*)V^(e) = 0, ^Ne{^PN,XN}. (4.22) 

where we used eq. (|4.3|) squared. This is just eq. (|4.9|) with a = 0, replacing p(0) by the universal 
parameter p(A*). From eqs. (|4.1U|) and ()4.11|) we see that the solutions are Bessel functions of index 
ib^, which are just sine and cosine. The differential equation no longer distinguishes between even 
and odd N, but the parity fixes the solution uniquely. The solutions for the asymptotic Cauchy 
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transform will no longer be singular as we expand around the point A = A*. However, we can still 
use the results of the previous subsection for a = in order to fix the solution uniquely depending 
on the sign 9m(e) in the Cauchy transform. Omitting the constant proportionality factor w{X^,)^^ 
when switching from wavefunctions to polynomials and their Cauchy transforms yields the explicit 
expressions for the latter: 

P2niO ~ cos(7r/)(A*)^) , 
P2n+i(.0 ~ sin(7r/5(A*)0 , (4.23) 



/i2n(C) 



exp(+i7rp(A*)C) , 9m(C) > 
- exp(— i7rp(A*)C) , 9m(C) < 



, / -i exp(+i7rp(A*)C) , 9m(C) > , . 

\-ieM-iMK)0 , 9m(C)<0 • (^-^^^ 

Due to the fact that no a-dependence is present in eq. ()4.22|) the polynomials and Cauchy transform 

in the chiral ensemble are identical to the even half of the non-chiral ones. The eqs. ()3.25|) and 1)3. 31() 
then give: 

Pn^kif) - cos(vr/,(A.)C) , (4.25) 
(a) ,^2^ / C"^exp(+i7rp(A*)C) , 9=m(C) > 



'^n,ch^^) ~ I _^"iexp(-i7rp(A,)C) , 9m(C)<0 ' ^ ' 

4.3 The soft edge 

Finally we turn our attention to the correlation functions in the vicinity of the edges A = ito of 
the support of the spectrum, the so-called soft edges [SUES]- Restricting ourselves to A = +a for 
symmetry reasons we define the microscopic scaling limit at the soft edge as 

ivi(A-a) = ^ . (4.27) 

In this way we magnify the region around the endpoint of the support. For negative (positive) ^ we 
are inside (outside) the support. In the limit eq. (|4.27|1 only the first term in eq. (|4.7|) contributes 
and the different power in is due to the factor (a^ — A^) that multiplies A{X)'^. In the large- limit 
the differential equation acquires the form: 

Vn{0-^^^Vn{0 = , ifNe {V-TV, Xn} . (4.28) 



It is again universal, extending the results of jl7j to the Cauchy transforms. After changing variables 
the solutions is obtained in terms of the Airy functions. Again omitting the prefactors w{a)^'2 when 
changing from wavefunctions to the original polynomials and Cauchy transforms we obtain: 

Pn{0 ~ K\{{2a-^A{afY i) , (4.29) 

/lAr(C) ~ < 



e+^Ai ((20^^(0)2)1(6+^) , 9m(C)>0 

(4.30) 



e- — Al[{2a-^A{af)-iQ(i-—] , 9m(C) < 
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Here we have chosen for the Cauchy transforms one of the two solutions of the Airy equation, that 
one which is hnearly independent of Ai(x) in coordination with the sign of 9m(^): 

Ai(^ze^^^ = ^e^f (Ai(z) =FiBi(z)) (4.31) 

The functions Bi(z) is a new feature emerging when negative moments of the characteristic poly- 
nomials are considered. For positive moments only (SHI, as well as for the spectral densities and related 
objects ^linHEl] only functions Ai(z) feature in the final results. 

The discontinuity between the two solutions is again a simple Airy function, due to the functional 
relation [55] : 

Ai{z) = -e+~Aifze+~j - e"~Aifze"~j (4.32) 

For real negative argument the Airy functions Ai and Bi are oscillating while for positive argument 
they monotonously decrease/increase, respectively. For large positive argument we conclude that 

lim =Fe=^^Ai fze^^^) = -iBi(z), (4.33) 

and thus see that both solutions eq. ()4.3U|) agree. This fixes our choice of the boundary condition as 
for large positive arguments outside the support the Cauchy transforms should no longer experience 
a discontinuity across the real line. 

The solution for even and odd is eqs. ()4.29() and 1)4. 3U() are the same. This immediately 
implies that the results for the UE and the chUE are identical, as must be the case far from the 
origin of the spectrum. Therefore we have to compute also the difference between Lp^q and (/^at-i as 
they both appear in expressions relating the correlation functions containing ratios of characteristic 
polynomials to orthogonal polynomials and their Cauchy transform'^. For that aim we use explicitly 
the differentiation eqs. 1)3. 8|) and 1)3. 20() for the corresponding wave functions together with the limiting 
relation eq. (j4.4|) . Following jSUl we finally obtain: 

</'Af-i(0 = <*^Af(6 H r^^ H^'nH) > "PN e {^AT, Xn} ■ (4.34) 

Ni.A{a) 

This completes our demonstration of universality for the orthogonal polynomials and their Cauchy 
transforms for both the UE and chUE, in all three different regions of the spectrum. 

We close with the following comment. Once the universality of expectation values of ratios of 
characteristic polynomials eq. (|2.2|1 is demonstrated it immediately translates to that of all other 
correlation functions that are generated by them by differentiations. As a check we can consider the 
correlation function of the eigenvalue densities by using eqs. ()4.13|) - ()4.14|) eqs. ()4.18|1 - 1)4. 19(1 respec- 
tively, and taking necessary derivatives and discontinuities. In this way one reproduces the known 
universal microscopic spectral densities of the unitary and chiral unitary ensemble, respectively, see 
e.g. ^Ij. It covers the bulk limit by setting a = 0. The difference is that in the chUE the resolvent has 
a true singularity at the origin while in the UE the discontinuity originates from patching two different 
solutions, each of which is regular at this point. This is another way of seeing that for the UE the 
"old" replica limit using saddle point techniques is only able to reproduce the asymptotic expansions 
|36j while for the chiral case it yields the exact result [SZj . Similar checks for the universal microscopic 
eigenvalue density at the soft edge can be performed using eqs. 1)4. 29() and ()4.3U|) . recovering known 
expression from ilTj . 

■^The same remark in fact applies already to the chUE in the previous subsections. 
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5 Discussion 



For a broad class of measures we have shown that the Cauchy transforms of orthogonal polynomials 
display properties very similar to those of the polynomials themselves. In particular, they obey 
the same three step recursion relation. The wave functions associated with the Cauchy transforms 
satisfy the same second order differential equation as those associated with the polynomials. This 
correspondence holds for finite-A^ with polynomial potentials of arbitrary but fixed even degree d, 
provided that N > d — 2. In the microscopic large- limit we have thus obtained the same universal 
Bessel equation at the hard edge, and the Airy equation at the soft edge of the spectrum. For the 
Cauchy transform the choice of the correct solutions is dictated by their singularity at the origin and 
the sign of the imaginary part of the argument. This identification was also supported by a comparison 
with exact integral representations available for the Gaussian weight at finite- A^. 

Using the general results IQi the universality of arbitrary ratios of characteristic polynomials 
in RMT with unitary invariance thus follows. Our consideration complements the rigorous proof 
|25j valid for the bulk of the spectrum and relied on the Riemann-Hilbert approach. The arguments 
presented here rely on the previous results in |3Ull21j . Although partly heuristic by strict mathematical 
standards, the method has the clear advantage of being much more transparent and able to treat both 
cases of the chUE and the UE on equal footing. 

Since our differential equation for the wavefunctions is valid at finite- A^ we could be tempted to take 
the macroscopic large- A^ limit as well. It gives the eigenvalue correlations at distances large compared 
to the mean level spacing, where the local fiuctuations are smoothed. While the spectral density itself 
is non-universal, all higher order connected correlation functions are known to be universal (Mj- It has 
been shown in that the density and the connected two point function can be also obtained from 
smoothing products of orthogonal polynomials. However, for higher order correlations the corrections 
to the orthogonal polynomials become extremely difficult to evaluate and are unknown so far. For 
that reason we have not tried to consider smoothened Cauchy transforms in order to recover the 
known one and two-point function. It is also less clear if the procedures of smoothening and taking 
the discontinuity after differentiation can be interchanged safely. 

Finally it would be very interesting to extend the general results |1U! to matrix models with complex 
eigenvalues and to prove their universality. First results in this direction have already been obtained 
j39j . There, arbitrary products of characteristic polynomials have been expressed in a determinant 
formula, which allowed to demonstrate their universality at the regime of weak non-Hermiticity. 
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and comments on the manuscript. This work was supported by Brunei University Vice Chancellor 
Grant (YVF) and by the Heisenberg fellowship of the Deutsche Forschungsgemeinschaft (GA). 



References 

[1] E.V. Shuryak and J.J.M. Verbaarschot, Nu cl Phys . A560 (1993) 306 |hep-th/9212088] ; J. Ver- 
baarschot, Phys. Rev. Lett. 72 (1994) 2531 |hep-t h79401059 . 

[2] A.V. Andreev and B.D. Simons, Phys. Rev. Lett. 75 (1995) 2304. 

[3] J. P. Keating and N. Snaith, Commun. Math. Phys. 214 (2000) 57. 

[4] E. Brezin and S. Hikami, Commun. Math. Phys. 214 (2000) 111 !math-ph/9910005| . 



16 



P.H. Damgaard, Phys. Lett. B424 (1998) 322 hep-th/9711110 ; G. Akemann and P.H. Damgaard, 
Nucl. Phys. B528 (1998) 411 hep-th/98 01133| ; Phys. Lett. B432 (1998) 390 ,hep-th/9802174, . 

J.C. Osborn, D. Toublan and J.J.M. Verbaarschot, Nucl. Phys. B540 (1999) 317 
[hep-th/9806110 ; P.H. Dam gaard, J.C. Osborn, D. Toublan and J.J.M. Verbaarschot, Nucl. 
Phys. B547 (1999) 305 hep-t h/9811212 . 

K. Splittorff and J.J.M. Verbaarschot Phys. Rev. Lett. 90, 041601 (2003) |cond-mat/0209594| . 

E. Kanzieper, Phys. Rev. Lett. 89, 250201 (2002) |cond-mat/0207745| . 

Y.V. Fyodorov and E. Strahov, Nucl. Phys. B647 (2002) 581 |hep-th/0205'215| . 

Y.V. Fyodorov and E. Strahov, math-ph/020405 l1. 

Y.V. Fyodorov and G. Akemann, (to appear in JETP Letters) |cond-mat/0210647| . 
M.L. Mehta, Random Matrices^ Second Edition, Academic Press, London 1991. 
E. Brezin and A. Zee, Nucl. Phys. B402 (1993) 613. 

G. Akemann, P.H. Damgaard, U. Magnea and S. Nishigaki, Nucl. Phys. B487 (1997) 721 
|hep-th/9609174|. 

E. Kanzieper and V. Freihkher, Phil. Mag. B77 (1998) 1161 1cond-mat/9704149| . 

P.H. Damgaard and S.M. Nishigaki, Nucl. Phys. B518 (1998) 495 |hep-th/971l'023l ; Phys. Rev. 
D57 (1998) 5299 hep-th/9711096 . 



E. Kanzieper and V. Freihkher, Phys. Rev. Lett. 78 (1997) 3806 chao-dyn/9701006| . 

G. A kemann, P.H. Damgaard, U. Magnea and S.M. Nishigaki, Nucl. Phys. B519 (1998) 628 



[hep^th/9712006| 

G. Akemann and G. Vernizzi, Nucl. Phys. B631 (2002) 471 |hep-th/020Tl65| . 
L. Pastur and M. Shcherbina, J. Stat. Phys. 86 (1997) 109. 

P. Deift, T. Kriecherbauer, K.T.-R. McLaughhn, S. Venakides and X. Zhou, Comm. Pure and 
Appl. Math. Vol. 52 (1999) 1335; ibid. 52 (1999) 1491. 

P. Bleher and A. Its, Ann. of Math. 150 (1999) 185 |math-ph/9907025l . 

P. Deift, Orthogonal Polynomials and Random Matrices (Courant Lecture Notes v. 3), Providence 
R.I., American Mathematical Society 2000. 

P.H. Damgaard, Proceedings of "Quantum Chromodynamics 98", Paris 1998, p. 296-303 
hep-th/9807026 . 



E. Strahov and Y.V. Fyodorov, |math-ph/0210010| 

K.B. Efetov, Supersymmetry in Disorder and Chaos, Cambridge University Press, Cambridge 
1997. 

[27] T. Guhr, J. Math. Phys. 32 (1991) 336. 



17 



[28] Y.V. Fyodorov, Nucl. Phys. B621 [PM] (2002) 643 |math-ph /0106006|. 

[29] P.J. Forrester, Nucl. Phys. B402 (1993) 709; T. Nagao and P.J. Forrester, Nucl. Phys. B435 
(1995) 401. 

[30] E. Kanzieper and V. Freilikher, NATO ASI, Series C (Math, and Phys. Sciences), Vol. 531, p. 
165 - 211, Kluwer, Dordrecht 1999 cond-mat/9809365j . 



[31] P. Bleher and B. Eynard, |^p-th/0209087| 

[32] C.A. Tracy and H. Widom, Commun. Math. Phys. 159 (1994) 151; ibid 161 (1994) 289; ibid 
163 (1994) 33. 



[33] E. Brezin and S. Hikami, Phys. Rev. E62 (2000) 3558 |math-ph/0004018. . 

[34] P.J. Forrester and N.S. Witte, Commun. Math. Phys. 219 (2001) 357 |math-ph/0103025| . 

[35] M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions, Dover Publ., New York 
1972. 

[36] D. Dalmazi and J.J.M Verbaarschot Nucl. Phys. B592 (2001) 419 |hep-th /0005"229l. 

[37] G. Akemann, D. Dalmazi, P.H. Damgaard and J.J.M. Verbaarschot, Nucl. Phys. B601 (2001) 77 
|hep-th/0011072| . 

[38] J. Ambj0rn, J. Jurkiewicz and Y. Makeenko, Phys. Lett. 251 (1990) 517; J. Ambj0rn, L. Chekhov, 
C. F. Kristjansen and Yu. Makeenko, Nucl. Phys. 404 (1993) 127 |hep-th/ 9302'oT4l. 

[39] G. Akemann and G. Vernizzi, Nucl. Phys. B in press |hep-th/021205l| . 



18 



